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We determine the energy-dependent amplitude for elastic tttt P-wave scattering in isospin-1 by 
computing part of the discrete energy spectrum of QCD in finite cubic boxes. We observe a rapidly 
rising phase shift that can be well described by a single p resonance. The spectrum is obtained 
from hadron correlators computed using lattice QCD with light quark masses corresponding to 
rriTT ~ 400 MeV. Variational analyses are performed with large bases of hadron interpolating fields 
including, as well as fermion bilinears that resemble qq constructions, also operators that resemble 
pairs of pions with definite relative and total momentum. We compute the spectrum for a range of 
center-of-mass momenta and in various irreducible representations of the relevant symmetry group. 
Hence we determine more than thirty values of the isospin-1 P-wave scattering phase shift in the 
elastic region, mapping out its energy dependence in unprecedented detail. 



PACS numbers: 14.40.Be, 12.38.Gc, 13.75.Lb 

I. INTRODUCTION 

Hadron spectroscopy is principally the study of reso- 
nances which decay strongly, with widths of tens or hun- 
dreds of MeV, into asymptotic states corresponding to 
a multiplicity of hadrons stable under the strong inter- 
action like the pion. Hadron resonances typically appear 
as enhancements in the continuous energy distribution of 
these multi-hadron final states. The simplest examples 
are elastic (or nearly elastic) resonances like the p(770) 
which appears in the / = 1, J = 1 (isospin-1, spin-1) 
channel of tttt scattering, or the A(1232) in the / = |, 
J = I channel of ttN scattering. In the elastic case the 
energy dependence of a partial- wave (definite- J) ampli- 
tude; c;an be expressed in terms of a single real number, 
the phase-shift, which in the case of a narrow resonance 
will show a rapid rise from angles near 0°, through 90°, 
approaching 180° at energies above the resonance. 

Although our only observables involve asymptotic 
many-hadron states, our desire is to understand hadron 
resonances at the level of interacting quarks and gluons 
within QCD, and this presents a significant theoretical 
challenge. The stable hadrons (e.g. pions) in the ini- 
tial and final states are strongly interacting quark-gluon 
composites and the interactions that give rise to the res- 
onant intermediate state contain nontrivial features like 
quark-antiquark annihilation. These complex and inter- 
related effects demand a consistent treatment of the non- 
perturbative dynamics of QCD. 

Lattice QCD, in which the theory is formulated on a 
cubic Euclidean space-time grid of finite extent, has the 
advantage of being a first-principles approach to QCD in 
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which the approximations required to render the theory 
computationally tractable are under control. It is a non- 
perturbative approach that, in the formulation we will 
use, respects unitarity - a particularly important feature 
since we plan to study resonances (like the p) which max- 
imally saturate unitarity. The lattice grid spacing acts 
as a regulator for the field theory which can be progres- 
sively reduced toward the continuum. In practice the use 
of 'improved' discretised actions leads to discretisation 
errors which are small for many quantities of interest. 
The restriction to finite spatial volume turns out to be 
a tool that allows us indirect access to hadron scattering 
amplitudes. These amplitudes are not directly accessible 
via (n > 2)-point Euclidean correlation functions [1] but 
can be inferred using the discrete energy spectrum which 
follows from periodic boundary conditions applied to a 
finite volume. The formalism for relativistic elastic scat- 
tering in a cubic box is presented in [2, 3] for the case of 
a system in its rest-frame, with the extension to moving 
frames in [4-6]. 

A well known practical problem with the implemen- 
tation of lattice QCD is the poor scaling of the required 
computation time for a realistic calculation with decreas- 
ing value of the quark mass. Only very recently have we 
seen calculations with quark masses low enough that the 
determined pion mass comes out near the experimental 
value, and they arc typically of only the simplest of quan- 
tities (see a review in [7]). In this paper we will present 
calculations performed with a single strange quark of ap- 
proximately the correct physical mass and two degen- 
erate light quarks with a mass such that the pion has 
m - 400 MeV. 

The challenge then is to compute the excited state 
spectrum of QCD in a finite-volume so that it can be re- 
lated to scattering amplitudes. In a series of papers [8-14] 
we have explored the problem of extracting excited state 
energy spectra from large matrices of 2-point correlation 
functions, (0|C'/(t)C'|(0)|0), where O are composite op- 
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orators with the quantum numbers of hadrons, built from 
quark and gluon fields. Using a large basis of such op- 
erators, wo have extracted highly excited spectra using 
variational analysis. Computing with such large operator 
bases is made efficient using the distillation framework 
[15], which also renders simple the inclusion of quark 
annihilation diagrams with a high degree of statistical 
precision. This was demonstrated in the computation of 
the isospin-0 meson spectrum in [13], where for example, 
the rj and rj' masses were determined with statistical un- 
certainty below 2% and even exotic isoscalar states with 
J^^ = 1 above 2100 McV could be cleanly extracted. 

In [8, 9, 11, 13], the meson operator basis was limited to 
constructions of the type V'f'V') where T was constructed 
using Dirac gamma matrices and up to three gauge- 
covariant derivatives. The resulting spectra did not show 
the expected strong volume dependence of multi-hadron- 
like states, and the proposed explanation was that the 
fermion bilinear operators have only weak overlap onto 
such states. To remedy this we should augment the op- 
erator basis with some having large overlap onto multi- 
hadron-like states. 

In [16, 17] we considered the case of the empirically 
non-resonant / = 2 tttt scattering using a basis of tttt- 
like interpolating fields of the form (i/)r^i/')pi ' ('0r7r'0)p2, 
the use of distillation allowing for operators of definite 
relative momentum at both source and sink and thus en- 
abling a variational analysis of matrices of correlation 
functions. By considering multiple frames in which the 
entire tttt system is in-flight and many irreducible repre- 
sentations of the reduced symmetry group of a boosted 
cube, the S and D-wave scattering phase-shifts were 
mapped out across the elastic region 2m7r < E^m < ^itIt^, 
showing the expected weak repulsive interaction. 

In this paper we will present results combining these 
two operator bases, forming matrices of correlation func- 
tions that include both constructions of the form ipTtp 
and (V^r^V)?! ■ (V'r'7rV')p2 with vector quantum numbers. 
The resulting variationally obtained spectrum leads to a 
detailed mapping of the P-wave phase shift for isospin-1 
TTTT scattering. 

The / = 1, J = 1 scattering channel has been con- 
sidered previously in lattice QCD calculations, [18- 22]. 
For example, in [19], four quark masses {m^ = 290 — 
480 MeV) were considered in two-flavor QCD calcula- 
tions, with up to six points on the phase-shift curve deter- 
mined at each pion mass. Fits to a simple Brcit-Wigner 
form were used to suggest a quark mass dependence 
for the p mass and width. In [20], a single pion mass, 
= 266 MeV in two-flavor QCD on a relatively small 
volume {rrinL ^ 2.7) was considered, with five points on 
the phase-shift curve determined; this was again fit with 
a simple Breit-Wigner form. One of our objectives in 
this paper is to map out the energy dependence through 
determination of many discrete values of the phase shift, 
conclusively demonstrating resonant behaviour and justi- 
fying a resonant parameterisation. This will be achieved 
through extraction of the spectrum on three volumes and 
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TABLE I. The lattice ensembles and propagators used in this 
paper. The lattice sizes and number of configurations are 

listed, as well as the number of time-sources (which varies 
somewhat according to the correlator momentum and irrep) 
and the number of distillation vectors TVvccs (to be described 
in Section VI) featuring in the correlator construction. 



for a large number of irreducible representations in mov- 
ing frames^. 

In each irreducible representation, we will extract a 
spectrum of states from threshold up to high excita- 
tions with many states above the inelastic thresholds into 
KK^TTUJ, Since we do not in this first study use in- 
terpolating fields that have good overlap onto these ad- 
ditional miilti-hadron channels, wo do not trust the de- 
termined spectrum outside the elastic region and restrict 
ourselves to the extraction of the elastic P-wave phase- 
shift between the tttt and KK thresholds. 

The remainder of the paper is organised as follows. In 
Section II we review the details of our lattice setup and in 
Section III we describe the method by which excited spec- 
tra are extracted from correlation functions. Section IV 
highlights the importance of including multi-hadron-like 
operators, in Section V we discuss the consequences of 
performing calculations on a lattice in a finite volume 
which has a reduced symmetry, and in Section VI wc 
summarise how correlators are constructed using the dis- 
tillation technique. We show the extracted finite-volume 
spectra in Section VII, in Section VIII wc present the 
energy-dependent phase shift which follows from these 
spectra, and we give a summary in Section IX. 



II. FINITE VOLUME LATTICE GAUGE FIELDS 

In Euclidean time, excited state contributions to cor- 
relation functions decay faster than the ground state, 
and at large times they are swamped by the signals of 
lower states, thus complicating the extraction of excited 
states. To ameliorate this problem we have adopted 
a dynamical anisotropic lattice formulation of Clover 
fermions with two light quarks and one strange quark. 
In this anisotropic formulation, the temporal extent is 
discretized with a finer lattice spacing than the spatial di- 
rections [26, 27], allowing a more precise resolution of the 
discrete time-dependence of correlation functions. In this 
work, computations were performed with spatial lattice 
spacing ~ 0.12 fm, and a temporal lattice spacing ap- 
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0.06906(13) 


atrriK 
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atirin 
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0.15678(41) 




3.444(6) 



TABLE II. Stable meson masses and the anisotropy ^ = as /at 
determined on the lattice ensembles listed in Table I. Pion 
and kaon masses come from an infinite-volume extrapolation, 
while the rj and u) masses are those evaluated on the 24^ lat- 
tice. 

proximately 3.5 times smaller, corresponding to a tempo- 
ral scale a^^ ^ 5.6 GcV. Results are presented for quark 
mass parameters a^m; = —0.0840 and ajTOs = —0.0743 
corresponding to a pion mass of 391 MeV, and on lattice 
sizes of 16^ X 128, 20^ x 128 and 24^ x 128 with correspond- 
ing spatial extents L ~ 2fm, ^ 2.5 fm and ^ 3fm. Some 
details of the lattices and propagators used for correla- 
tion constructions are provided in Table I. The masses 
of the lightest stable mesons containing light and strange 
quarks, as determined on these lattices, are shown in Ta- 
ble II. 

As well as the volume dependence of energies of multi- 
hadron states originating in hadronic interactions which 
are the subject of this work, there can also be depen- 
dence of single- hadron energies on L. These polarization 
effects, corresponding to a single hadron 'sensing' itself 
around the periodic volume, are typically largest for the 
lightest hadron, the pion, and can be characterized by 
the product ttIt^L. At the single quark mass value used 
in this study, m,ri' ranges from 3.8 to 5.7, and since ef- 
fects typically fall off exponentially with increasing m^L 
we expect them to be rather small here. As an example, 
previous investigations have found [16, 28] that the pion 
mass variation with m^j^L on these lattices is not large. 

In an infinite volume of continuous, isotropic space- 
time, a single free particle can have any value of momen- 
tum with magnitude varying continuously up from zero, 
with the energy related to the momentum by the rela- 
tivistic dispersion relation E'i = + In a finite 

L X L X L volume with periodic spatial boundary condi- 
tions, the momentum of a free particle is restricted to the 
discrete set p= ^ [ux, Uy, Uz) = ^n, with the n's being 
a triplet of integers. We will write momenta in units of 
^ with square brackets, e.g. p = [nj;,ny,nz]- 

A further complication which arises from our use of an 
anisotropic lattice is the need to determine the precise 
value of the anisotropy, ^, which relates the spatial lattice 
spacing to the temporal lattice spacing at = a^/^. The 
anisotropy appears in the dispersion relation of a free- 
particle, where as found previously [17], the dispersion 
relation for the pion (and other stable hadrons) can be 
well described by a continuum-like form 

{a^Enr = {am^)' + ^(^^y\nf, (1) 



used throughout the rest of this paper. 

III. EXTRACTING EXCITED STATE SPECTRA 

Our determination of the spectrum of eigenstates of 
QCD in a finite volume proceeds through the calcula- 
tion of matrices of correlation functions between suitable 
hadronic creation and annihilation operators at time 

and t respectively, 

a,(t) = (o|o,(t)o](o)|o), 

For each correlation matrix, the set of operators {Oj}, 
constructed from color-singlet combinations of quark and 
gluon fields, all have the same conserved quantum num- 
bers. Within the basis of operators used we can at- 
tempt to find the optimal linear combination for inter- 
polation of each possible finite- volume eigcnstatc |n). A 
commonly-used method to achieve this is a variational 
solution [2, 29, 30] and herein our particular applica- 
tion of the variational method follows that developed in 
Refs. [8, 9, 31]. A system of generalized eigenvalue equa- 
tions is established for the correlation matrix 

Cit)v^{t) = Xr^{t)C{to)v^{t), (2) 

where and v" are eigenvalues and eigenvectors for a 
state labelled by n. Eq. (2) is solved for the eigenvalues, 
and the exponential dependence on the Euclidean time, 
An(i) ^ p--En(t-*o)^ used to determine the energy En 
of the state. The orthogonal^ eigenvectors represent the 
optimal combination of the operators Oi for interpolation 

of the state |n) from the vacuum, (On^* )^ ^ '^fC'l- 

Any two-point correlation function on a finite spatial 
lattice can be expressed as a spectral decomposition 

n " 

where the "spectral overlap factors", 
Z" = (n|C||0) arc related to the eigenvectors by 

Z'; = y/kSnC^"^"/^ v'J*Cji{to). The form of Eq. 3 as- 
sumes that t ^ T, the temporal length of the box, so 
that the contributions arising from other time orderings 
on the periodic lattice can be ignored. However, in 
practice, these alternate time-ordered contributions have 
a small but discernible effect on the determined energies; 
a method to account for them was presented in Ref. [17] 
and is used again in this work. 

IV. INCLUDING MULTI-HADRON 
OPERATORS 

Although we have previously reported on the extrac- 
tion of a significant number of excited meson states 
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FIG. 1. The extracted spectrum of states, on the 24^ volume, using variational analysis on a matrix of correlation functions with 
total momentum P — [0, 0, 1] transforming irreducibly in the Ai representation of the relevant group. The three columns left to 
right correspond to a large basis featuring operators of both fermion bilinear structure and tttt structure, just fermion bilinear 
and just tttt respectively. The histograms show the relative size of matrix element overlaps (nlO^IO) for a range of operators. The 
red, gray and blue bars are fermion bilinears while the orange bars (drawn somewhat wider to be more visible) are tttt operators. 

Also shown by the dashed orange lines are the energies of non-interacting tttt pairs, E^'^' — \J + nf ( x ) ^ + \l + ^2 
and the dashed purple line shows the position of the KK threshold above which the system is inelastic. 



[8, 9, 11, 13, 14], those calculations used only operators 
built from a single fermion bilinear at rest or in flight, 
Oi ~ "ipTiip. A large basis was built using up to three 
gauge-covariant derivatives acting on the fermion fields, 
and this allowed many states to be extracted from varia- 
tional analysis of correlator matrices. It was found that 
the extracted spectrum, which could be largely described 
in terms of constituent qq constructions supplemented 
with 'hybrid' states in which a gluonic excitation is also 
present, showed only very weak dependence on the lattice 
volume. In particular, the spectrum lacked signatures of 
multi-meson states, which if non-interacting would ap- 
pear with a strongly volume-dependent spectrum and a 
characteristic distribution across irreducible representa- 
tions of the lattice symmetry group. In the interacting 
theory we expect the actual eigenstates to be (volume- 
dependent) admixtures of what we might call "single- 
hadron" basis states and multi-hadron basis states. 



Our explanation for this observation of a restricted 
spectrum is that the fermion bilinear operators have an 
overlap onto multi-hadron states that is suppressed by 
powers of the lattice volume, such that they couple only 
very weakly into the multi-hadron sector. To the ex- 
tent that the eigenstates in any particular volume are 
admixtures of single-hadron and multi-hadron states, one 
can argue that correlators computed using fermion bilin- 
ears contain contributions from all eigenstates, |n), of 

the appropriate symmetry in ^2e' e"-^"*. This is 
true, but there remains the practical problem of extract- 
ing the spectrum, which is particularly challenging when 
there are nearly-degenerate states, a likely occurrence in 
a dense spectrum of hadron scattering states. An effec- 
tive solution, which does not depend upon distinguishing 
small energy differences, is to make use of orthogonality 
within the variational method, as described in the pre- 
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FIG. 2. EfTective mass of Eq. 4 with atEi = 0.1654, atE2 — 0.1779 for a range of mixing angles, 8. Left panel shows typical 
time-separations in our calculation over which the time-dependence is essentially flat. Right panel indicates that eventually the 
presence of two low-lying states could be detected but only through observation of unrealistically long time separations with 
high statistical precision. 
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FIG. 3. Effective mass of lowest principal correlator obtained 
from variational analysis using only fermion bilinear opera- 
tors (more time sources were used here than presented in the 
central column of Fig. 1). The gray band indicates a fit to 
the principal correlator using two exponentials as described in 
the text. The dashed gray line is the lower energy in that fit. 
Red and purple bands indicate the energies obtained from the 
lowest two principal correlators when using the full operator 
basis (with a smaller number of time sources). 



vious section, where each state is optimally interpolated 
by a different orthogonal linear superposition of the basis 
operators even if the states are approximately degener- 
ate. 

The lack of operators that have significant overlap on 
multi-hadron states causes the variational analysis to 
be unable to form orthogonal eigenvectors in the larger 
space {lhlS?!r„)'lhadron)}' sincc they have access only to 
{Ihadron)}- prcscnt a simplc model below of what can 
happen if only a restricted operator basis is used when 
the true eigenstates of the system in finite- volume are ad- 



mixtures of single-hadron and multi-hadron basis states. 

In [17] we demonstrated the construction of TTTr-like 
operators built from products of variationally optimised 
pion operators of definite momentum, TT^{k). In simple 
terms we construct 

fcl ,fe2 

for definite total momentum P — ki + k2 in irreducible 
representation, A. The sum is over the different direc- 
tions, for fixed magnitude, of the definite pion momenta 
fci,fc2, with the weights C ensuring that the operators 
transform in A. Through the use of multiple possible 
choices of magnitudes, |fci|, |fc2|i we can build a basis of 
operators. 

The importance of including TTTr-like operators in ad- 
dition to fermion bilinears is demonstrated in Fig. 1 
which shows the extracted spectrum computed on the 
24^ lattice when the system has a total momentum of 
P — [0, 0, 1] and lies in the irreducible representation Ai, 
roughly corresponding to helicity-0 (to be discussed more 
in the next section). The three columns show the spec- 
trum extracted using the full basis of operators (fermion 
bilinears plus TTTr-like), only fermion bilinear operators, 
and only TTTr-like operators. There are striking differ- 
ences, notably the large gap between atE^^ = 0.20 and 
0.30 which, when tttt operators are included, is populated 
by several states which have large overlap onto the tttt 
operators. For our current purpose, the most important 
difference is the presence of two states in the elastic re- 
gion with the full basis while only one appears with the 
restricted operator bases. Notice that in the full basis 
these two states each have large overlap onto both the 
fermion bilinear operators and the tttt operator of lowest 
relative momentum. 

The essential problem with using a restricted opera- 
tor basis can be demonstrated using a simple two-state 
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mixing hypothesis, where the state basis is a "would-be" 
stable single hadron state, jp), and a multi-hadron state, 
Itttt). The finite- volume energy eigenstates are linear su- 
perpositions of these two basis states: 

\Ei) = cose\p) + sin6'|7r7r) 
\E2) = - sin 6*1/9) +cos6'|7r7r). 

A variational analysis that used multiple operators with 
good overlap onto both \p) and Itttt) would be able to re- 
solve the two orthogonal combinations, with two separate 
principal correlators having time dependences ^ e~^i*, 
~ e~^^*. However, if one restricts to only operators with 
good overlap onto |p), the variational method will not 
be able to form two orthogonal vectors and the principal 
correlator will have a timc-dcpcndcnce similar to 

(p|e-^*|p) =cos2 0e-^i* + sin2 6'e-^^*, (4) 

which features both energy scales. In practice these en- 
ergies can be very close together, e.g. atEi = 0.1654(7) 
and atE2 = 0.1779(7) giving a* AS = 0.013(6) in the case 
presented in Fig. 1. In Fig. 2 we show the effective mass^ 
corresponding to Eq. 4 for a range of mixing angles, 9. 
One clearly sees that over the time-separations typically 
resolved in lattice QCD calculations and, considering typ- 
ical statistical uncertainty, the time-variation is compati- 
ble with a single-state hypotheses. It is only at very large 
times that the correlator relaxes to the true ground state 
and the correlator can be seen to clearly contain more 
than one state at a low energy scale. In practice correla- 
tors will also contain excited-state pollution from higher 
energy scales that impact at small Euclidean times. 

In Fig. 3 we show the effective mass of the lowest 
principal correlator obtained from variational analysis us- 
ing only the formion bilinear operators. The gray curve 
shows the effective mass corresponding to a fit to the 
principal correlator, between t/at = 5 and t/at = 41, of 
the form X{t) = (1 - A)e"-^(*~*«) + Ae^-^'^*"*") where 
atE = 0.16793(28) and UtE' = 0.494(11). The xV^dof 
for the fit is below 1.0 indicating that the data is de- 
scribed well in terms of one low-lying state and ex- 
cited state contributions at a much higher energy scale. 
Clearly, from this principal correlator alone, we cannot 
infer the presence of two low-lying energy levels at 0.1654 
and 0.1779. 

In summary, if one wishes to determine the spectrum 
reliably, it proves necessary to include explicit multi- 
hadron operators into a variational basis. In this paper 
we will consider the inclusion of operators resembling a 
pair of pions of definite total and relative momentum. 
Since we do not yet include any operators resembling 
pairs of kaons, the next lightest stable hadrons, we will 
be cautious about trusting extracted levels near the top 
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of the elastic region where they may be mixing with KK 
basis states. 

In the next section we present in somewhat more detail 
the construction of two-meson operators relevant for a 
study of TTTT scattering in isospin-1. 



V. MESON-MESON OPERATORS ON A FINITE 
CUBIC LATTICE 

The symmetry of a lattice with a finite extent is re- 
duced compared to that of an infinite volume contin- 
uum. In our implementation, we have a cubic lattice 
discretisation in a cubic box with periodic boundary con- 
ditions. The appropriate symmetry group is the dou- 
ble cover of the octahedral, or cubic, group with parity, 
O^, and this is the symmetry relevant for a system of 
hadrons overall at rest. The allowed momenta arc quan- 
tised by the boundary conditions giving p= ^ (n, m, p) , 
OT p = [n,m,p] in our notation, where L is the spatial 
extent of the lattice (in physical units) and n,m,p arc 
integers. For a system 'in flight', i.e. with a non-zero 
overall momentum, P ^ 0, the symmetry is reduced fur- 
ther to that of the little group [32], which wc denote by 
LG(P), the subgroup of under which P is invariant. 

The consequences of this reduced symmetry were dis- 
cussed in detail in Ref. [17] and here we briefly review the 
salient points relevant for tttt scattering in isospin-1. For 
P = 0, the continuum spin, J, is no longer a good quan- 
tum number on the lattice and states arc instead labelled 
by the irreducible representations, irreps, of the octahe- 
dral group, of which there are a finite number. Parity, 
P, and any relevant flavour quantum numbers are still 
good. The manner in which the various components of 
a state are distributed, or subduced, into the irreps, 
A^, is presented in Table II of Ref. [17]. 

The situation is more complicated for P 7^ 0; the pat- 
tern of subductions of the various helicities, A, into the 
little group irreps depends on LG(P), i.e. the type of 
momentum P, see Ref. [11]. Table II of Ref. [17] shows 
these subductions. Note that, apart from the A = com- 
ponents, parity is not a good quantum number for a sys- 
tem in flight, but any relevant flavour quantum numbers 
are still good. 

A state consisting of two identical hadrons must have 
a definite symmetry under the interchange of the two 
hadrons; it must be symmetric under this interchange 
if the hadrons are bosons or antisymmetric if they are 
fermions. This generalises to systems of two hadrons, 
having identical masses, related by a symmetry, e.g. Tr+, 
TT~ and 7T^ if isospin is a good symmetry, or K and K 
related by charge-conjugation. In particular, taking into 
account the symmetry of the spin and any flavour parts 
of the wavefunction, the spatial wavefunction must have 
some definite symmetry under this interchange and so is 
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P LG(P) 




[0,0,0] 0° 


Tf 1\ 3' 
Tj- 3I 
^2 


[0, 0, n] Dic4 


^1 1\ 3^ 
£2 1\ 32 
Bi 3I 
S2 3I 


rn T) Til OiTo 


yli 1\ 3^ 
A2 31 
Bi 1\ 32 
B2 1\ 32 


[n, n, n] Dies 


A2 3' 
£2 l\ 3^ 


[n, in, 0] „ 
[n, n, mj 


^1 1", 3* 
yl2 l\ 3^ 



TABLE III. The pattern of subductions of / = 1 tttt partial 
waves, < 3, into lattice irreps, A, where N is the number 
of embeddings of this i in this irrep. This table is derived 
from Table II of Ref. [17] by considering the subductions of 
the I when P = Q and the subductions of the various holicity 
components for each I when P 7^ 0. Here P is given in units 
of and n, m are non-zero integers with n ^ m. We show 
the double-cover groups but only give the irreps relevant for 
integer spin. 

restricted to either only even or only odd partial waves*, 
For TTTT in isospin-1, the flavour (isospin) wavefunction 
is antisymmetric and so Bose symmetry restricts the sys- 
tem to odd partial waves (with parity P = (—1)^ = — 1 
when P = 0). The reduced symmetry of the finite vol- 
ume lattice means that multiple partial waves appear in 
a single irrep as shown in Table III. 

We note that, in contrast, for a two-meson state such as 
tttt', where tt' represents an excited pion, or ttw, there is 
in general no such Bose symmetry constraint. For such 
a system overall at rest, parity, P oc (—1)^, is a good 
quantum number and so, even on a finite volume lattice, 
odd and even partial waves do not appear together in any 
irrep"''. However, when P 7^ 0, parity is not a good quan- 
tum number, there is in general no definite symmetry 
under the interchange of the two hadrons^ and all (. can 
mix; there is no separation between odd and even partial 
waves. In the current work, we restrict ourselves to scat- 
tering below inelastic thresholds and so do not consider 



* e.g. if isospin is a good symmetry, even (odd) tttt partial waves 
have isospin / = or 2 (i = 1) with positive (negative) G-parity. 
KK states with positive (negative) eharge-eonjugation parity, C, 
~ [KK + C KK), only occur in even (odd) partial waves 

® e.g. non-interacting 7r(fc)cj(— fc) and 7r(— fc)aj(fc) have the same 
energy and appropriate linear combinations arc eigenstates of 
parity 

® e.g. non-intera<;ting 7r(fci)a)(fc2) and 7r(fc2)w(fei) with P = k\-\-k2 
do not in general have the same energy 



such mixing between odd and even partial waves. This 
will however be relevant as we go higher up in energy and 
in other scattering channels. 

We use the single and multi-meson operators, respect- 
ing the symmetries of the finite volume lattice, which 
have been developed in a series of papers [8, 9, 11, 17]. 
For a single meson at rest, we first construct an oper- 
ator with definite continuum and J^-component M, 
Qj ,M(^j^ = 0)^ consisting of a fermion bilinear featuring 
gauge-covariant derivatives and Dirac gamma matrices, 
and with the correct flavour structure. We then sub- 
duce these into the relevant octahedral group irreps to 
form lattice operators [8, 9]. At non-zero momentum, 
our construction considers the subduction of an opera- 
tor with definite continuum helicity, A, O"' '^(fc), where 
k is quantised as discussed above, into the relevant little 
group irreps [11]. The general pattern of subductions is 
summarised in Table II of Ref. [17]. 

In Ref. [17] we discussed the single-meson operators 
relevant for a pion, which has 7=1 and negative G- 
parity, at rest and in flight. At rest, the subduction of a 
= 0~ operator into the irrep A^ = is trivial, 

C>[°:1(0) = O0"(0) . 

^1 

At non-zero momentum, for all the momenta we consider, 
the subduction into the A2 irrep is also trivial. 

When we use the variational method to flnd the optimal 
linear combination of operators to interpolate a pion, we 
will include operators subduced into from other J 
(for P — 0) and other helicity. A, subduced into A2 (for 
P 7^ 0). We use ir^k) as a shorthand to represent these 
optimal operators in the appropriate irreps. 

For the p meson, which has 1=1 and positive G- 
parity, at rest the three-dimensional Tj~ irrep is a simple 
representation of = 1~ and the subduction is straight- 
forward; this is given explicitly for our basis choice in 
Appendix A of Ref. [9]. In flight, the A = components 
subduce into the Ai irrep and, depending on LG(P), 
the A = ±1 components subduce into either the two- 
dimensional E2 irrep or the one-dimensional Bi and B2 
irreps [11]. 

As mentioned above and discussed in detail in Ref. [11], 
for a system in flight, parity is in general not a good 
quantum number. The hclicity-0 components have a 're- 
flection parity', rj = P(— 1)"', where J and P are respec- 
tively the spin and parity at rest. For the momenta we 
consider, A = with 77 = -|-(— ) subduces into the ^41(^2) 
irrep and so there is no mixing between these even on a 
flnite volume lattice. However, for A 7^ there is no such 
symmetry and therefore, for example, the A = ±1 com- 
ponents of the p{l~) appear in the same lattice irreps as 
the A = ±1 components of the 61 (1^) (both these mesons 
have the same flavour structure). This is not relevant in 
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the current work where we are restricting ourselves to rel- 
atively low energies, below inelastic thresholds, but must 
be taken into account if we wish to consider higher en- 
ergies. Similar situations occur in other scattering chan- 
nels. 

In Section IV we have argued that to reliably deter- 
mine the complete spectrum we need operators that have 
good overlap onto multi-hadron basis states. We find 
that operators built from the product of two or more 
single-hadron operators of the form described above are 
able to achieve this. We will follow the construction used 
in Ref. [17] where a general tttt creation operator is, 



, Jfei 
(tttt) - . 



fc2]t 



fcie{fei}* 
fee {__&}* 



(5) 

where 7r(fc) is a single-pion operator and C is a Clebsch- 
Gordan coefficient for Ai (g) A2 — )• A with Ai 2 = of 
O)^ if fci,2 = and Ai,2 = A2 of LG(fci,2) if ki,2 0, and 
where A is an irrep of LG(P). The sum over fci.2 is a sum 
over all momenta in the stars of ki^2, which wc denote by 
{^1,2}*, and by which we mean all momenta related to 
ki,2 by an allowed lattice rotation. The ki, k2 label on 
the TTTT operator indicates that it was constructed from 
single-pion operators with momenta in {ki}* and {^2}*- 
We refer to Ref. [17] for further details, a discussion of 
the Clebsch-Gordan coefficients, and explicit values of C. 

As discussed in Ref. [17], we use "optimised" single- 
pion operators, 7r(fc), in order to obtain tttt correlators 
which are dominated by ground state pions at smaller 
times, i.e. to reduce contamination from tttt', where tt' 
is any higher mass state with pion quantum numbers. 
We follow that reference and construct variationally op- 
timal operators by diagonalising a matrix of correlators 
in large basis of operators for each relevant irrep. As 
described in Section III, the n"^ eigenvector gives the 
variationally optimal linear combination of the basis op- 
erators to overlap with state n. Our basis consists of 
all single-meson operators subducing into the relevant 
irreps constructed from any possible Dirac gamma ma- 
trix structure and 0,1,2,3 derivatives at rest [9] or 0,1,2 
derivatives in flight [11]. The efficacy of these operators 
was demonstrated in Ref. [17] and their use allows us to 
perform analyses at smaller Euclidean times. 

In the current work we restrict ourselves to overall mo- 
menta P = [0,0,0], [0,0,1], [0,1,1], [1,1,1] and [0,0,2]. 
The various combinations of fci and ^2 used in our tttt op- 
erator constructions are presented in Table IV. For each 
lattice irrep, we include in our basis these tttt operators 
along with all relevant isospin-1 single-meson operators 
from any possible Dirac gamma matrix structure com- 
bined with 0,1,2,3 derivatives at rest [9] or 0,1,2 deriva- 
tives in flight [11]. The number of operators we use in 
each irrep is summarised in Table V. 



p 


Volumes 


ki 


k2 


A(^) 


[0,0,0] 

^ tl 


16^,20^24=' 


[0,0,1] 
[0, 1, 1] 


[0,0,-1] 
[0,-1,-1] 


Tr 


16^ 


[1,1,1] 


[-1,-1,-1] 


Tr 


[0,0, 1] 
Dic4 


16^20^24^ 


[0,0,0] 
[0,-1,0] 

[-1,-1,0] 
[0,0,-1] 


[0,0, 1] 
[0, 1, 1] 

[1,1,1] 

[0,0,2] 


Ai 

Ai,E2, Bi * 
A\,E2,B2 * 
^1 


[0,1,1] 
Dic2 


16^,20^24^ 


[0,0,0] 
[0,1,0] 
[-1,0,0] 
[1,1,0] 
\0 1 -11 


[0, 1, 1] 
[0,0,1] 

[1,1,1] 

[-1,0,1] 
21 


Ai 
Bi 

Ar,B2 
Bi, B2 
A-i Si 


[1,1,1] 

Dies 


16^20^24^ 


[0,0,0] 
[1,0,0] 


[1,1,1] 

[0, 1, 1] 


^1 

Ai,E2 




[2,0,0] 


[-1,1,1] 


Ai,E2 


[0,0,2] 


16^,20^24^ 


[0,0,0] 


[0,0,2] 


Ai 



TABLE IV. The two-pion operators used presented for each P 
on various volumes; also shown is LG(P). Example momenta 
fei and k2 are shown; all momenta in {ki}* and {fe}* are 
summed over in Eq. 5. Swapping around ki and gives the 
same operators up to an overall phase. (*these Bi and B2 
were not considered on the 16^ volume) 



P Irrep 


Single-meson tttt 20^ 24* (16^) 


[0,0,0] T- 


26 2 (3) 


Ai 

[0,0,1] g 
B2 


18 4 (4) 
29 2 (2) 
9 1 (0) 

9 1 (0) 


Ai 

[0, 1, 1] Bi 
B2 


27 3 (3) 
29 3 (3) 
29 2 (2) 




21 3 (3) 
35 2 (2) 


[0,0,2] Ai 


18 1 (1) 



TABLE V. Number of single-meson and tttt operators used 
for each P and irrep on the various volumes. The two-meson 

operators are listed in Table IV and all relevant single-meson 
operator structures are considered including up to 3 deriva- 
tives at rest and up to 2 derivatives at non-zero momentum. 
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VI. CORRELATOR CONSTRUCTION 
THROUGH DISTILLATION 

In Section IV we have emphasised the need to include 
operators that have a strong overlap with expected multi- 
meson states and the desire to perform a variational anal- 
ysis of a matrix of correlation functions forces us to find 
a correlator construction method that allows for such op- 
erators at both source and sink. In the isospin-1 chan- 
nel, the contractions to form these correlators will in- 
volve quark-line 'annihilation' - that is diagrams which 
feature propagation of a quark from t to the same t. In 
the previous section we described how multi-hadron op- 
erators transforming irreducibly under the relevant sym- 
metry can be constructed from products of single-hadron 
operators of definite momentum. The projection into def- 
inite relative momentum requires a sampling of all spatial 
sites on a timeslice. The correlator construction method 
known as distillation [15] satisfies all the above desider- 
ata in a natural and efficient way. 

Distillation is a quark-field smearing method which is 
designed to increase overlap onto the low modes relevant 
for low-lying hadronic states - we define a smearing op- 
erator on a time-slice which acts in 3-space (x) and color 

(*) 

N 

□ (B, yj; t)^Yl ^) ^n(yj; t), (6) 

n=l 

where one option is to choose {^„} to be the lowest N 
eigenvectors of the gauge-covariant Laplacian on times- 
lice t. If all the quark fields in a correlator are smeared by 
application of this operator, the combination of eigenvec- 
tors and the Dirac matrix inverse, M"^, called a 'peram- 
bulator', SX{t')M~'^{t' ,t)S,mit) = Tnm{t',t) wiU appear 
in any propagation. Thus the basic numerical problem 
to be solved is inversion of the Dirac matrix on sources, 
{^n}n=i...N, which is a smaller vector space than that of 
the full lattice. A detailed presentation of the properties 
of distillation can be found in Ref. [15]. 

Five basic topologies of diagram appear when one Wick 
contracts correlators in a basis of ■ipT-ip and ^pT^^ ■ {pT^tp 
operators projected into overall isospin-1 - these are 
shown in Fig. 4. All five topologies require light-quark 
perambulators between fixed igrc and varying t, and these 
are obtained by inversion of the Dirac matrix from a 
source on tsrc- While the propagation from tgrc to tgrc 
seen in topologies ^42 and A44 does not require any ex- 
tra inversions, the propagation from t to t seen in A24 
and A44 requires inversion from sources at every t we 
wish to consider. Such T{t,t) perambulators, computed 
for both light and strange quarks on the 16^ lattice, were 
previously used in [13] in the computation of isoscalar 
correlators featuring rpTip operators. In this paper we 
use T(t,t) perambulators computed on 16'^, 20^, 24'^ lat- 
tices for all 128 timeslices. 

In order to reduce statistical noise, and make 
more use of each configuration, we compute correla- 




FIG. 4. A schematic of the basic Wick-contraction topologies 
required to compute a matrix of correlation functions in a 
{V>ri/), %l)T^%l) ■ tpT'^Tp] basis with isospin-1. tsrc is fixed and 
the time dependence with respect to t — tsrc is studied. 

tors for several values of tsrc and average over them, 
C{t) = Sig,^ C{t-\- tsrc, tsrc)- Although we have in- 
verted from every timeslice on the lattice (in order 
to have the T{t,t) perambulators), covariance between 
neighbouring timeslices reduces the effectiveness of aver- 
aging over isrc once the sources get too close together, 
so in practice, given the computational cost of contract- 
ing perambulators into correlators, we average only over 
a handful of well-separated tsrc values, typically fewer 
than 8 per configuration. 

As an example of the distillation contraction for 
^pT^tp ■ -tjjT^ip at both source and sink, choosing the 
temporal origin to be at tgrc = 0, we have 

where the 3-space (x) and color (j) dependence of 
Tit)xi,yj might include gauge-covariant spatial deriva- 
tives and/or a projection into definite momentum. In- 
tegration of the quark fields leads to terms which include 
the Wick contraction ^144, 

A44 = r„„(o,i)<i>;^„(t)T„p(i,0$f,(t) 

XT,,(i,0)C(0)T.„(0,0)$f„(0), 

where the repeated distillation indices are summed over 
and where the Dirac-spin indices are suppressed. This 
factorised form is such that any operator construction, 
r (embedded in = ^t^£.m), can be considered. For 
example, if there is a momentum projection, a Dirac spin 
projection (F"^) and gauge-covariant derivatives {D), 
then 

xi,yj 
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5 10 15 20 25 30 t/aj 



FIG. 5. P = [001], Ai correlators evaluated on the 24^ 
lattice constructed from tp"fitp and the lowest momentum 
TVTv operator ~ 7r([001])7r([000]). Shown is e'^^C{t) with 

atE = atEthr. = \j (2afm,r)'^ + (^' ^Lja ) • Normalised to all 

be equal to 1 at t/at = 31. The inset shows the data at 
t/at = 31 indicating the relative statistical noise on the cor- 
relators. 



and provided e*^'^ TD . . .1 ^. ^.{t) is well-supported in the 

restricted space of vectors {£.n{xi)]^_^ ^, distillation 
should provide an efficient method to compute the corre- 
lator. In practice [9, 11, 15, 17] we find that the modest 
numbers of vectors, N, given in Table I, are sufficient 
to obtain excellent signals for hadrons with momenta 

We demonstrate the quality of the determined cor- 
relators in Figs. 5 and 6. In Fig. 5 we present corre- 
lators, evaluated on the 24'^ lattice in the Ai irrep of 
P — [001], constructed using the simplest fermion bilin- 
ear operator, tl^^iil^, and the lowest momentum tttt opera- 
tor 7r([001])7r([000]). We see that correlators featuring 
multi-pion operators are not much noisier than simple 
fermion bilinears and that both types of operator are 
likely to have a significant overlap with the ground state 
and differing overlaps with excited states. In Fig. 6 we 
show that the statistical noise on correlators does not 
grow rapidly with increase in the momentum of pions in 
TTTT operator constructions. 



VII. FINITE- VOLUME SPECTRA 

We computed correlation matrices for a range of to- 
tal momentum P and irreps A using the fermion-bilinear 
plus TTTT operator basis described in Section V. The corre- 
sponding finite- volume spectra were determined by appli- 
cation of the variational method described in Section III. 
As an illustration, we show in Fig. 7 the 'principal cor- 
relators', An(t), for the lowest eight levels (from a total 



FIG. 6. P = [001], Ai correlators evaluated on the 24^ lat- 
tice. The statistical noise, relative to signal, is shown for tttt 
diagonal correlators for a set of pion momenta. 

of 22) in the case P = [001], A = Ai on the 24^ lattice. 
The fits to the time-dependences of these determine the 
energies with high precision. The corresponding opera- 
tor overlap structure was previously displayed in the first 
column of Fig. 1. The effective masses of the five lowest 
principal correlators are shown together in the bottom- 
right pane of Fig. 7. 

In Fig. 8 we show the extracted volume dependence 
of the spectra for P = [001] for the irreps, Ai and E2- 
Also shown are the energy thresholds for various inelas- 
tic processes and the energies of non-interacting meson- 
meson states TTTT, KK, ttuj. Clearly we are not observing 
the expected number of levels in the inelastic region and 
we expect that this is caused by our not using operator 
constructions with good overlap onto e.g. KK. As men- 
tioned in Section V, in the E2 irrep, parity is not a good 
quantum number and contributions from — 1+ can 
appear - indeed the points shown in grey at atE^m ^ 0.25 
appear to have large overlap with operators characteris- 
tic of the bi meson. Since the operator basis we have 
used does not sample well the inelastic spectrum, we will 
restrict our analysis to the elastic region between tttt and 
KK thresholds. 

In Fig. 9 we show the volume dependence of the spec- 
tra across P and A in the elastic region. For the phase- 
shift analysis which follows in the next section we ex- 
clude points very close to KK threshold since we have 
not included KK-like operators in the basis. Levels 
determined on the 16'^ lattice, which potentially suffer 
from the largest unwanted exponential mjrL dependence 
and which are often noisy, are mostly excluded from the 
phase-shift analysis. 
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FIG. 7. Result of variational analysis of P = [001], Ai correlation matrix. For the lowest eight levels, we show the "principal 
correlators", the eigenvalues \n{t) of the variational method, Eq. 2. Plotted as e^"'*~*°' An(f) are the data and the fit according 
to the form An(i) = (1 — An) e~^" + An e~^"^*~*°K Energies are those measured in the rest frame of the lattice. 
The bottom-right pane shows effective masses of X„{t) for the lowest five levels. Horizontal lines indicate the fitted energies, 
at En 




FIG. 8. Volume dependence of P = [001] spectrum in irreps Ai and E2. Black and gray points show the extracted spectrum. 
Dashed lines show energy thresholds and solid lines the non-interacting meson-pair energies. 
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FIG. 9. Volume dependence of elastic spectra for various P, A. Plotted is atE^m versus L/ag. Also shown by dashed horizontal 
lines are the tttt and KK energy thresholds. Solid curves indicate the non-interacting tttt energy levels. Points shown in gray 
are excluded from the phase-shift analysis in the following section. 
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VIII. ELASTIC SCATTERING PHASE SHIFT 

Our implementation of the formalism relating elastic 

TTTT scattering phase-shifts to the finite-volume spectrum 
at rest and in flight is presented in Section VIII of [17]. 
Here we repeat only the essential feature, that the phase- 
shifts are related to the finite-volume energies by a for- 
mula which for tttt in / = 1 takes the form, 



det 









„2i53(-Ecm) 



-u 



P,A 



LV 



= 0. 

(7) 



The matrices here, in the space of £ contributing to the 
irrep, arc formally infinite dimensional, with U being a 
matrix of known functions particular to the little group 
irrep (-P, A). Which values of £ contribute to a given 
(P, A) is determined by the subductions given in Ta- 
ble III. We presume we can truncate the matrices to a 
finite size since at low energy we expect 5£>3 ^ ^ (5i . 

In Fig. 10 we present the determined phase-shifts as- 
suming that all de>s are negligible throughout the elastic 



region. In this case Eq. 7 becomes one equation in one un- 
known, di{Ecm), for each energy level, Ef.m{P, A, n) in the 
elastic region. The assumption 5^>i will be shown 
to be justified in Subsection VIII B. 



A. Resonant Parameterisations 

The phase-shift in Fig. 10 shows clear resonant be- 
havior in a small region around aji^cm = 0.152 which 
suggests we might attempt to describe 5i{Ecm) using a 
paramcterisation featuring a single resonance. A popular 
choice is a variant of the relativistic Breit-Wigner with an 
energy-dependent width enforcing the p^m^^ behavior as 
Pern mandated by angular momentum conservation: 



tan(5i(i;cm) 



Ecm ri(i?cm) 



m% 



where T^^^{E^n.) 



67r 



(8) 



The parameterisation of the width is in terms of a cou- 
pling constant, g, which, it has been suggested, is largely 
independent of the quark mass [19, 33]. 



We choose to fit the energy levels directly within a function of the form, 
X\{ai}) = J2 J2 [^cm(i; PAn) - i?P- (i; P^n; {aj)] C-\L; PAn; P'A'n') [E,ru{L; P'A'n') - £;p- (L; P'AV; {aj) 

■f- PAn 
P'A'n' 

(9) 

where E^^^- {L; PAn; {ai}) is the n**^ solution of Eq. 7 with a parameterised phase-shift depending upon parameters 

with values {flj}, e.g. the resonance mass m,R and the coupling, g. Data covariancc, C, whose off-diagonal elements 
between energies evaluated on the same ensemble can be non-zero, is estimated using jackknife. 

I 



Fitting energy levels with a^i^cm < 0.18 to the Breit- 
Wigner form in Eq. 8, we obtain the following parameter 
values, errors and correlation. 



atmR 
9 



: 0.15241(33)(10) 
:4.83(13)(2) 



1 



-0.09 
1 



29-2 



1.81, 



where the first error is statistical and the second error 
refiects variation of m^r and ^ within their respective un- 
certainties (see [17]). There is observed to be very little 
correlation between the mass and the coupling constant. 
The energy dependence of the resulting 5i(i?cm) is shown 
by the blue curve in Fig. 11. 

A criticism of the parameterisation in Eq. 8 is that 
the form, introduced to give the right threshold be- 
haviour, continues to grow at large energies in an unreal- 
istic way. As seen in Fig. 11, the phase-shift approaches 
180° very slowly, allowing the effect of the resonance to be 
felt many half-widths above the mass. Generally one ex- 
pects the behaviour to be damped at larger energies. 



and one physically motivated approach is to appeal to 
the idea that the resonance has a finite spatial size, and 
introduce 'barrier factors'. One implementation, com- 
monly used in experimental partial wave analysis, follows 
von Hippcl & Quigg [34], by considering the interaction 
giving rise to the resonance to have a sharp size, R. The 
resulting scattering wavefunctions are outgoing spherical 
Hankel functions which give rise to polynomial damping 
factors. In the case £ = 1 this leads to the following 
modification of the energy-dependent width: 



g' pI^ 1 + {pnRf 



(10) 



where pn is the cm scattering momentum at the reso- 
nance mass, PR = \\/m\ — 4m^. The barrier factor is 
normalised such that it does not affect the width eval- 
uated at the resonance mass. We fit the same data as 
above to obtain 
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FIG. 10. P-wave tttt elastic scattering phase-shift, 5\{Ecm), determined from solution of Eq. 7 applied to the finite-volume 
spectra shown in Fig. 9 under the assumption that (5^>i = 0. Energy region plotted is from tttt threshold to KK threshold. 



atruR = 0.15226(34)(11) 

g = 5.06(15)(2) 
R/at = 16.6(52)(17) 



-0.14 
1 



-0.09 
0.32 
1 



dof 



43,6 



= 1.68, 



which shows a sUghtly improved quality of fit, although 
there is clearly some correlation between the coupling g 
and the range R. The range expressed in physical units 
R = ^^ii^ 0.6 ± 0.2 fm would seem to be reason- 



able on the usual hadronic scale. The resulting energy 
dependence is shown by the red curve in Fig. 11 where 
it is seen to approach 180° more rapidly than the simple 
Breit-Wigner. 

The particular form of the damping function is a 
model-dependent choice and we can explore the sensi- 
tivity by trying other parameterisations. For example a 
gaussian form (previously considered in a quark model 
study [35]), 



9 Pc 



(11) 



Fitting the same dataset we obtain 

atmji = 0.15224(34)(14) 
g = 5.08(17)(3) 
at(3 = 0.029(7)(3) 



1 -0.18 0.16 
1 -0.47 
1 



dof 



43,5 



1.67, 



indicating that the particular functional form of the 
damping appears to be relatively unimportant. In phys- 
ical units, /3 = at/3 ■ ^^^^ w 160(40) MeV. The energy 



dependence is shown by an orange curve in Fig. 11 that 
lies almost exactly on the red curve already described. 

Another parameterisation that has been used to fit ex- 
perimental phase-shift data is provided by Pelaez and 
Yndurain (see Ref. [36] and their subsequent papers). 



cot 5i{Ea 



2p; 



2ml 
mj^Ec, 



Ei 



Bo + B: 



Ecm - y/gp - E^^ 



which, while it appears cosmetically to be very differ- 
ent to a Breit-Wigner, in fact has an energy depen- 
dence which is rather similar, with the three parameters 
•niR, Ba, Bi able to conspire to provide damping. The 
additional parameter, sq, is not allowed to float, and fol- 
lowing the proposers' suggestion is set to 2-01^^ + Wp, as 
determined on this lattice, atJ~SQ = 0.29. Fitting yields 



atruR = 0.15227(34)(12) 
Bo = 2.71(77)(21) 
Bi = 6.0(33)(9) 



1 



-0.06 
1 



-0.05 
0.99 
1 



dof 



43,7 



1.68, 



a reasonable description of the data. The extremely high 
degree of correlation between Bq and Bi suggests that 
they may not be the most natural way to parameterise 
this amplitude. The energy dependence is plotted in 
Fig. 11 using a green curve that lies almost exactly under 
the orange and red curves already plotted. 

We have presented the data and fits in units of 
the inverse temporal lattice spacing thus far to avoid 
ambiguities with how one sets the lattice scale. If 
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FIG. 11. P-wave tttt elastic scattering phase-shift, <5i(i5cm), as determined by describing the finite-volume spectra by resonant 
parameterisations as described in the text. The barrier factor variations and the Pelaez & Yndurain fit all lie on top of each 
other. Also shown as gray points are the data previously presented in Fig. 10. Energy region plotted is from tttt threshold to 
KK threshold. 



we choose our usual scale setting procedure where 
o-t = usine; the Vl baryon mass determined 

on these lattices {atm^ = 0.2951) and the physical fl 
baryon mass rriQ^^^ — 1672 MeV, then the simple Breit- 
Wigner fit corresponds to m^; = 863.5(19) (6) MeV and 
Tfi = 10.1(6)(1) MeV. As expected in a calculation with 
heavier than physical mass light quarks, the resonance 
mass is somewhat larger than the physical p mass. The 
small width is explained by the much-reduced phase- 
space for decay of an 864 MeV resonance into two pions 
of mass 391 MeV compared with the physical kinemat- 
ics. We observe from the rather similar x^/^dof that 
the data do not clearly distinguish between the various 
parameterisations which vary only in the tails of the res- 
onance. This may be due to the very narrow nature of 
the resonance with the small phase space for decay. 



B. Role of higher partial- waves 

From Eq. 7 and Table III it is apparent that in principle 
many partial waves contribute to the determination of 
the finite- volume spectrum in each irrep, in particular 
when the system is in-fiight. The next lowest i that can 
contribute in tttt J = 1 scattering is £ — 3 which is leading 
in e.g. irreps (P = [001], Bi) and (P = [001], B2). For 
the lattice volumes we consider, the lowest energy level 
in these irreps is always above the elastic region, and 
as such we cannot apply Eq. 7 without concern about 
neglecting other open channels (in this case KK). If we 
assume that there is zero coupling into KK and proceed 



in a cavalier manner with application of Eq. 7 we obtain 
points at energies only slightly above the KK threshold 
that have 63 compatible with zero (roughly (—1 ± 1)°). 

One way to obtain estimates of S3 in the elastic 
regime is to consider a number of approximately de- 
generate energy levels coming from different irreps. By 
writing a version of Eq. 7 for each one we can ap- 
proximately solve that coupled set of equations for 
Si , ^3 at the relevant energy. This approach was de- 
scribed in some detail for tttt / = 2 scattering in 
[17]. An example set of levels is ([000], Tf, n = 0), 
([001],i;2,n = 0) and ([Oil], B2, n = 0) which on the 24^ 
lattice all have an energy OtE^cm ~ 0.153(1). Solving the 
coupled system of equations we find Si = 145.7(22)° and 
S3 = -0.048(55)°. The same set of levels on the 20^ lat- 
tice have atEcm ~ 0.155(1) and give Si = 151.1(30)° and 
S3 = -1-0.002(24)°. 

We also tried parameterised fits to all data points, as 
in the previous section but including a scattering length 
parameterisation for the £ = 3 wave, pl^ cot ^3 = 1/03, as 
well as a resonant parameterisation of <5i(ii^cm)- The fits 
were of essentially the same quality (in l^doi) and gave 
03 — —3.4(33) (6) X 10^ -aj with a negligible change in the 
£ = \ Breit-Wigner parameters. This parameterisation 
gives ^3 = -1.3(13)° at the KK threshold. 

In summary, the lattice data require no non-zero value 
of ^3 throughout the elastic region and our analysis in 
the previous section based upon 53 = is justified. 

These observations (at ^ 400 MeV) are in ac- 
cord with experimental expectations (at the physical pion 
mass). In the tttt partial wave analysis of Estabrooks 
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and Martin [37], no F-wave amplitude was required 
in the elastic region to describe ■K~p — > Tr^Tr+n data. 
Above inelastic threshold, the i^-wave amplitude between 
Ecm = 1.5 and 1.9 GeV was well described by a Breit- 
Wigner form with barrier factors. Subsequent experi- 
ments have determined the resonance parameters of the 
/93(1690) with greater precision [38]. The tail of this reso- 
nance in the elastic region (which must vanish like p^m) is 
such that (^3 ^ +1.8° at the KK threshold. The energy- 
dependent phase-shift analysis of Tr+p — >■ tt+tt" A++ data 
by Protopoescu et al. [39] suggests a ^3 which is small and 
negative below KK threshold, with the largest deviation 
from zero ~ —1.5(5)°. 

IX. SUMMARY 

In summary, we have used lattice field theory methods 
to compute part of the discrete energy spectrum of QCD 
in finite boxcis. Through the known connection of the dis- 
crete spectrum to elastic scattering amplitudes, we have 
mapped out the energy dependence of the tttt isospin-1 
P-wave scattering phase shift in unprecedented detail up 
to the KK threshold. The data, presented in Fig. 10, un- 
ambiguously show the rapidly-rising form expected in the 
presence of an elastic resonance. As observed in Fig. 11, 
the phase shift can be well described by a single p reso- 
nance. While in principle the discrete spectrum can be 
sensitive to scattering in higher-partial waves, we find 
that no non-zero value of F-wave (or higher) phase-shift 
is required in the elastic region. 

The P-wave scattering phase-shift is summarised in 
Fig. 12, where the scale in all dimensionful quantities is 
set using ot = . A simple Breit-Wigner fit describes 

the data reasonably well in terms of a single narrow res- 
onance. We note that the extracted coupling g is com- 
patible with other lattice determinations [18-22]. 

This work represents a successful application in a res- 
onant channel of the methodology for lattice computa- 
tions of scattering phase shifts presented in Ref. [17], 
where it was initially applied to a non-resonant channel. 
Distillation enabled us to efficiently compute correlators 
with large bases of carefully constructed fermion-bilinear 
and TTTT-like interpolating operators, in various irreducible 
representations of the relevant symmetry group and for 
a range of center-of-mass momenta, with high statisti- 
cal precision. In addition, the inclusion of the necessary 
quark annihilation contributions is rendered straightfor- 



ward. Variational analyses of the resulting correlators 
gave a large number of finite- volume energy levels which, 
in turn, allowed us to determine an unparalleled number 
of points on the phase shift curve. 

We plan to compute scattering in other channels using 
the same technology, including situations where exper- 
imental and phenomenological understanding is incom- 
plete and where we cannot assume simple parameteri- 
sations will describe the data. In these cases, mapping 
out the phase shift in detail will be more important than 
in the simple case of elastic tttt scattering in isospin-1 
that we have considered here, where a single vector reso- 
nance is expected to dominate the scattering amplitude. 
Determining the energy dependence of the scattering am- 
plitude in detail will inform the parametcrisations. Fu- 
ture work to explore the hadron resonance spectrum will 
need to consider inelastic scattering - we can reliably ex- 
tract the finite-volume spectrum above inelastic thresh- 
olds by including the appropriate meson-meson-like op- 
erators (e.g. KK . . .), following the general methodology 
developed in Ref. [17]. 
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